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Abstract 

We propose an analysis of the dead time losses in counting imaging 
detectors such as MWPC which can resolve k simultaneous hits, and an- 
alyze in more detail an X — Y detector which has a third wire set which 
allows for the recognition of simultaneous impacts. 



1 Introduction: the standard X — Y detector. 

A particle detector with an intrinsic X — Y structure, such as a MWPC (see for 
example [3]) can be used as a pulse counting 2-D imaging detector. A specific 
pixel of the imaging detector shares a specific readout channel X\ with other 
pixels, and shares another specific readout channel y± with other pixels, but in 
such a way that the combination of X\ and y\ are not shared by any other pixel. 
A detector with n x x n y pixels can hence be read out with n x +n v channels. As 
long as there is a single hit on the detector, giving rise to a single X and a single 
Y channel being hit, the hit pixel can be reconstructed without ambiguity. As 
long as the signal development time of a single hit is small compared to the 
time between impacts, a time correlation of X and Y hits can associate the 
right A-hit with the right Y hit, and within this time window, a single hit will 
be able to be resolved into a single pixel. 

However, from the moment that two (or more) successive impacts on the 
detector cannot be resolved in time, we have two (or more) A-hits, say channels 
Xi and x%, and two (or more) K-hits, say channels y\ and y%, and it is not 
clear whether the association should be {x\,y\) and (2:2,2/2) or whether the 
association should be (£1,1/2) and (x2,yi) — if there is no other information, 
such as amplitude correlations, and the hits are similar. In other words, a double 
impact generates 4 possible pixel candidates, of which 2 are the correct ones, 
and 2 others are "ghost" hits. Let us consider that the largest time which can 
separate the A and the Y signal of a single hit, equals r, which we call the 
coincidence time. The only way to associate unambiguously an X hit to a Y hit 
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is that during the time r before and during the time r after the event, no other 
event occurs. The probability in a Poissonian stream with n hits per second on 
average, that during a time T no event occurs, equals 

P{T, k = 0) = exp (-nT) (1) 

The condition for an event, occurring at t — to be acceptable (no ghosts), is 
that during the interval t = — r to t = CP, no event happens, and that during 
the interval t = + to t = r, no event happens. These probabilities being 
independent for a Poisson stream (as the intervals are disjoint), the probability 
for this to be so equals P(r,k = 0) 2 = exp (— n2r). As such, the number of 
recorded hits (with no ambiguity) in an X — Y detector equals: 

m = nexp(— 2nr) (2) 

In other words, an X — Y detector behaves as a paralyzable (see for instance, 
[I]) detector, but with twice the coincidence time r as dead time. 

Of course, there is also an intrinsic dead time for each X wire itself, and 
each Y wire itself, which can often be longer than the coincidence time window 
r. If we call this wire dead time, t w , the time during which no two distinct 
hits on the same X wire can be distinguished, this will also give rise to a dead 
time, in parallel with the coincidence dead time above. In order to estimate its 
importance, let us suppose that our detector is uniformly irradiated. This means 
that an X wire has an average Poissonian flux of N x = —. For a given event at 
time t = 0, we don't want its X wire to be "blinded" by an event preceding it, 
outside of the already considered coincidence window. The probability that no 
event occurs on the same X wire in the time window starting a time t = —t w but 
stopping at t = — t, the start of the coincidence window, equals, using equation 
[U exp (—N x (t w — t)). We can apply the same reasoning for the Y wire, so the 
joint probability that no event occurs, nor on the hit X wire, nor on the hit Y 
wire, so that both the X and the Y wire are not "blinded" before the actual 
coincidence window opens, equals: 

P w = exp [ - — (t w - t) ) exp ( — —{t w - t) ) (3) 
\ n x J \ n y J 

The overall relationship between the true flux n, and the observed flux to, in 
the case of a uniform, Poissonian illumination, is: 



to = ncxp (— 2nr) exp ( (t w — r) ) exp ( (t w — r) I = rtexp {—nrt) 



n 



(4) 



where 



n = t I 2 - ) + t w [ — + — ) (5) 



n x n,, 



's, 

As such, a standard X — Y detector has a counting behaviour which is that of a 
paralyzable detector with dead time given by Tt. We also see that this dead time 
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consists of the "coincidence dead time" (2r) and has a contribution of the "wire 
dead time" of the order of 2r w /n x . So, the wire dead time is less important 
(under uniform irradiation) when we have r w /n x <C r, which is often the case 
for large MWPC. 

2 A third "disambiguation" electrode grid. 

An MWPC is often made by "sandwiching" an anode wire plane between two 
cathode wire planes, or between a cathode wire plane and a solid conductor 
surface. One can make use of this degree of freedom to use one cathode plane 
and the anode plane as the X — Y channels, and the second cathode as a third 
set of channels. In the case of a wire plane, one can imagine for instance a 
set of wires under 45 degrees with the X — Y grid. This has been proposed 
by Lewis [5]. It is also possible to use a more symmetrical setup where the 3 
wire planes (two cathode planes and one anode plane in the case of a MWPC, 
or 3 different readout directions on other types of detectors) make angles of 60 
degrees with each other, giving rise to some hexagonal 'honeycomb' structure. 
This idea has been implemented using multi-GEM detectors for X-ray imaging 
[6] , and the same idea has also been proposed and experimented for Chcrenkov 
photon detectors (where there are typically simultaneous hits), in [7]. In order 
to analyze the principle however, which we set out to do in this paper, the exact 
geometry doesn't matter, as long as the topological relationships between the 
coordinates are the same. 

Let us call this third set of channels, the Z channels. To each possible hit 
corresponds, as before, an Xi and a yj, but now also a z& channel hit, and 
the interesting point is that not all combinations of x, y and z correspond to 
existing pixels. In other words, for a single hit, there is some redundancy in 
the information. This redundancy can be used to try to find out what are the 
correct, and what are the "ghost" combinations of x and y, when we have a 
multiple-hit event. One can easily establish that it is always possible to find 
a numbering scheme of the channels X, Y and Z, such that, with a hit (x,y), 
there corresponds a Z-hit given by z — x + y. It is herein that lies the possibility 
to distinguish potentially multiple hits: of a list of X signals {xi,X2, ■■■}, a list 
of Y signals {2/1,2/2, ■■■} and a list of Z signals {zi, z%, ...}, not all combinations 
are possible and one can hope that only the correct combinations are allowed 
for. 

2.1 Two simultaneous hits. 

In fact, amongst two hits, the disambiguation is complete. Imagine the hits 
{2:1,052}) {2/1,2/2} and {21,22} (where we don't know of course the right num- 
bering when we reconstruct them). The right hits are (xi, 2/1, Z\) and (#2, 2/2, z-i)- 
This means that Z\ = Xi + yi and z 2 = x 2 +2/2- The question is: are there other 
possibilities ? Let us assume that ii ^ S2i !/i ^ 2/2 and z\ ^ z 2 for starters. 
Imagine that (xi, 2/2,^1) is a solution too. This means that Z\ — X\ +2/2, from 
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which follows that y\ — y^ what was against the starting hypothesis. A similar 
conclusion can be drawn for the case (x±, 7/2, z 2), and all other thinkable cases 
can be reduced to these two by permutations. On the other hand, imagine 
that x\ = X2- In that case, both the (xi,yi,zi) and (xi, 2/2,^2) are the right 
solutions, and there are no others (no ghosts can be constructed). 

So this means that two simultaneous hits can be recognized correctly with 
certainty. If we limit ourselves to this result, then we can calculate the re- 
lation between true counting rate and observed counting rate. An event will 
be lost, if there is more than one (other) event in the time frame that starts 
time r before and ends time r after it. The Poisson distribution is given by 

P(k, A) = cxp( ~, A)A " , from which it follows that P(k <u,X) = P(fc, A) = 

F p"„^ = Q(u, A) with Q the regularized incomplete gamma function (in [3J, the 
incomplete gamma function P(a,x) = 1 — Q(a,x) is introduced ; see also [4]). 
So the probability that an event gets accepted is equal to P(k < 2, 2nr) = 
r(2,2nr) = e~ 2nT (l + 2nr), which is the probability that in the interval of 
length 2t, centered onto the event under discussion (but without counting this 
event of course) there is one or no hits. As such, the observed counting rate will 
be: 

m = ner 2nT (1 + 2nr) (6) 

We see that (as far as coincidence dead time is concerned), a detector that can 
discriminate up to two simultaneous hits, does not exactly follow a paralysable 
or non-paralysable model, but does in fact, much better: there is no first order 
dead time! 

m ~ n - 2r 2 n 3 + ... (7) 



2.2 k simultaneous hits. 

Let us imagine for a moment that we have a detector that can discriminate, 
without any difficulty, up to k simultaneous events. What's the "dead time" 
now ? For a Poisson distribution, in order for an event to be counted, we need to 
have less than k events in an interval of 2t centered on an event to be counted. 
We will assume that the dead time is caused by the coincidence time window 
needed, and is not dominated by individual "wire" dead times. By a similar 
reasoning as above, we will obtain that the observed counting rate will be: 

fc_1 (2nT) u 

m = nP(< k, 2nr) = nQ(k, 2nr) = ne" 2 '"" ^ ^ f- (8) 

u=0 U ' 

For k = 2, we find the same expression as in equation [S] For k = 3, we find: 

m = ne- 2nT (l + 2nr + 2n 2 r 2 ) (9) 

This is interesting, because we see that the last term neutralizes the significant 
correction term in equation [7] We now have for low rates: 



771 — 11 — 



(10) 
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Figure 1: Observed counting rates as a function of incoming flux, when up to k 
events can be discriminated in a time r. 



In figure [TJ we show the relationship between a (normalized on 1/r) incoming 
Poissonian flux and the (also normalized) counting rate of discriminated events 
for k = 1,2,3,4,5. Often, the " dead-time correction" (which is given by 1—m/n) 
shouldn't be too important if one wants to give quantitative credibility to an 
imaging detector, and at the ILL, we take as a definition of acceptable counting 
rate for a detector, the one where the dead time correction is equal to 10%. The 
"dead time correction factor" 1 — m/n is shown in figure [5] In our case, we 
find then that this counting rate n w % is given by the solution of the following 
equation: 

Q(fc,2n 10% r) =0.9 (11) 
which leads to a solution using the inverse regularized gamma function (see [1]): 

n w % = y ( 12 ) 

Of course, the 10% can be judged a bit too severe, and some prefer allowing for 
20% correction, in which case we have: 

"20% = 7T ( 13 ) 

Solving these equations for the first 5 fc-values, we find: 
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Figure 2: Dead time correction factor as a function of incoming flux, when up 
to k events can be discriminated in a time r. 



k 


2n 10 %r 


2^20% T 


1 


0.105 


0.223 


2 


0.531 


0.824 


3 


1.102 


1.535 


4 


1.745 


2.297 


5 


2.433 


3.090 



So, for a traditional wire chamber, we find that a 10% dead time correction is 
reached when the flux is 0.105 times l/(2r); when the chamber can distinguish 2 
simultaneous hits, this flux is 5 times greater. When the chamber can distinguish 
3 simultaneous hits, this flux is 2 times greater again, etc... If we adhere to 
the definition of maximum flux at 20% dead time correction, being able to 
distinguish 2 simultaneous hits gives about 4 times higher counting rates and 
being able to distinguish 3 simultaneous hits adds a small factor 2 on top of 
this. 

We hence see that the biggest gain in counting rate (as of our definition 
using a maximum dead time correction of 10% or 20%) occurs when we go 
from single-hit to double-hit identification, which can, as shown previously, be 
obtained by using a third grid. 
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2.3 How many hits can a third grid accept simultaneously 

? 

If there are more than two simultaneous hits, then a third electrode cannot 
guarantee that no ghost hit will be present. Let us consider the case of 3 
simultaneous hits, and let us consider that (x\, 1/2, 23) is a ghost hit. If there is 
a ghost hit, we can always permute indices to have the X value be x\ and the Y 
value be 3/2- In this case, the Z value can only be Z3, as any other value would 
bring us back to the case of two hits and a ghost, which we demonstrated, is 
not possible. So the condition to have a ghost hit for a simultaneous impact of 
3 hits is: 

x\ + V2 = x 3 + y 3 (14) 

(or a similar condition for permuted indices). Nevertheless, most of the time 
when we have 3 hits the above condition will not be satisfied. In that case, the 
signals still allow for a non-ambiguous reconstruction of the 3 events. In order 
to find out in detail what percentage of possible fc-hit events can be resolved 
one can take two roads: one is, for a specific setup, to have a Monte-Carlo 
simulation of fc-hit events ; the other is to use some analytical estimations. In 
any case, the percentage of such resolved hits will depend also on the specific 
image that is projected onto the detector: we will assume uniform irradiation 
in our analysis. 

2.3.1 Analytical estimation of resolution of fc-hit events. 

Imagine that we have a fc-hit event. This means, a priori (we'll come to that), 
that we have fc X values, fc Y values and k Z values. We will assume (which 
is, depending on the exact geometry of the detector, only approximately true), 
that each of the n x X values are equally probable, that each of the n y Y values 
are equally probable, and that each of the n z Z values are equally probable. For 
each possible combination of an X value (in the hit list), and a Y value (also 
in the hit list), we can calculate the corresponding Z value (using z — x + y). 
If this Z value is in the list of Z hits, we have to accept the hit. Of course, for 
the real hits, this will be the case. It is for the fc(fc — 1) wrong combinations of 
X and Y values (xi,yj) that we shouldn't, by coincidence, fall on a Z value z m 
which is also present in the hit list. By a previous reasoning, we already know 
that this z m value cannot be the Z value which goes with the correct event of 
the X value Xi, or of the Y value yj, so the potential index m cannot be i or 
j. As such, there arc (fc — 2) possible Z values which could, by coincidence, 
be equal to Xi + yj. In order for our ghost couple (xi,yj) to be rejected, none 
of these should be equal to Xi + yj. Assuming (which is an approximation, 
but a reasonable one) that we can consider these (fc — 2) values (out of n z 
possible) as being statistically independent, the probability for Xi + yj to be 
equal to one of them equals then , so the probability not to have a ghost hit 
for (xi,yj) equals 1 — h=l. There are k(k — 1) of these different combinations 
(xi,yj) possible, so the probability (if, again, we consider all these potential 
coincidences as statistically independent, which is of course approximate but a 
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reasonable hypothesis) that none of these combinations gives rise to a ghost hit 
(and hence, the fc-hit event is totally identifiable), is then: 

^H'-^r) (15) 

In the above formula, we have put n w = n x = n y — n z the number of channels, 
chosen to be equal for X, Y and Z. What is clear is that when n w goes to infinity 
(infinitely many channels per coordinate), that all fc-fold hits are resolvable. It 
is only due to a finite number of channels that 3 or more hits can potentially 
give rise to ghost hits. 

There is a caveat, however. If we have fc hits and there is a finite number 
of channels, then there is also a finite probability that there are less than fc 
X- wires hit, or less than k Y wires hit or less than k Z wires hit because of the 
finite probability to hit the same wire twice. 

In the appendix in subsection IA.1[ the number Z l k is introduced, which 
gives us the number of different ways one can construct an ordered list of k 
elements, of which each element is one of n possible ones, and in which there 
are exactly I different elements present. This corresponds to the number of ways 
one can distribute k hits over n different wires, and touch in all I different wires. 
Under a uniform irradiation with "distinguishable" hits, each of these different 
ways is equally probable, and hence the probability, if there are A; hits, to have 
I wires hit, is given by: 

Z l 

PQ;k,n w ) = ^f^ (16) 

So a better approximation for the probability of being able to identify cor- 
rectly a fc-hit event, is to use formula [T5l with the average value of I: 

k 

{l) k = Y J P^k,n w )l (17) 
1=1 

to give: 

Pgood(k) =1 (18) 



2.3.2 Monte Carlo estimation of resolution of fc-hits 

Let us consider a detector with 3 wire grids, at 120 degrees one from the other 
and 32 wires each, equally spaced. We delimit the "useful" space as the hexag- 
onal surface that is covered by the 3 different wire sets simultaneously, which 
consists of 768 pixels (out of the 1024 crossing points of each pair of 32 x 32 wire 
planes). To each of these pixels correspond hence 3 numbers, %i,yi, Zi = Xi +yi- 
A fc-hit is formed by drawing independently fc natural numbers between 1 and 
768 (corresponding to the pixels), and looking up what x, y, z values they corre- 
spond to, to make up the lists of X, Y and Z hits. Next, all values of (xi+yj)^ 
arc tested against all values of z u with u these are the ghosts that cannot 
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be rejected ; we count how many there are. We repeat this procedure a large 
number of times (number of trials N) and make then a histogram of the number 
of ghosts that each trial generated. Zero ghosts means that the event could be 
correctly reconstructed. The fraction of the number of zero ghosts No(k) over 
N is an estimator of the probability to be able to reconstruct a fc-hit event. 
For N = 10000, we find: 
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If we compare that with our analytical estimation in equation 1181 (multiplied 
with N), then we obtain: 
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which gives quite good agreement especially for the lower k values. For the 
higher k values, there is an under-estimation of the number of correctly identified 
fc-hits. 

Repeating the experience with 128 wires per plane, we find, after 10000 trial 
events, for the Monte Carlo result: 
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while the analytical estimation gives us: 
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Again, one observes a relatively good performance for the lower fc-values, 
and an underestimation of the number of correctly identified events at higher k 
values. 



2.4 Count rates, dead times, and the third coordinate. 

Previously, we considered detectors which could, with 100% certainty, discrimi- 
nate an event when there were no more than k — 1 other events in a time slot 2r 
centered on our event, and which couldn't handle, also with certainty, an event 
when there were k or more events in the given time slot. The detector with a 
third coordinate however, has a different behavior: it can discriminate, with a 
certain probability P goo d{k), the case where there are (k — 1) events in the said 
time slot (on top of our event-under-test). This means that the relationship 
between incoming flux and observed, identified rate of hits is now given by: 

oo 

m = n 

P Poia(k - 1, 2 T n)P good {k) (19) 

fc=i 
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Figure 3: The observed (normalized) counting rate of identified hits as a function 
of the (normalized) incoming flux, using equation [T5] 

a sum which we can of course truncate to the first few terms given that as well 
the probability to have a fc-hit as well as the probability to resolve it correctly, 
will drop fast for high k values. As such, our analytical estimate can be used, 
given that its performance for relatively low k values is adequate. Using equation 
[T91 we can plot (figure [3]) the observed (normalized) counting rate of identified 
events as a function of the incoming flux, for different numbers of wires per 
plane. The associated dead time correction is shown in figure [4] Calculating 
the incoming fluxes that give rise to a dead time correction of 10% or 20%, we 
find (normalized onto 2r as was the case in the fixed k case): 
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2.734 
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3.467 
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Figure 4: The dead time correction of identified hits as a function of the (nor- 
malized) incoming flux, using equation 1191 



3 Conclusion. 

In this paper we discussed the theoretical dead time correction occurring in 
detectors which need a certain time to recognize the correct positioning of an 
impact in the case that they can handle k impacts simultaneously. We applied 
this to the specific case of a detector with an X — Y — Z planar wire structure. 
It is established that in principle, such a detector (which has 50% more readout 
channels than a standard X — Y detector) can accept incoming counting rates 
which are about 12 times (or more) larger than the standard X — Y detector at 
10% dead time correction. As such, the investment in 50% more channels can 
result potentially in 1200% higher counting rates. However, these results are 
only established in the case of a Poissonian, uniform irradiation. 

We established a relatively simple analytical estimate of the dead time cor- 
rection, and verified its applicability with a more detailed Monte Carlo simula- 
tion. 
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A Counting multi-hits. 



A.l Counting distinguishable hits. 

Let us define a fc-fold multiple hit event on a set of "pixels" with n possible 
"positions" as an ordered series {x\, X2, ...,Xk}, where Xi takes its values out of 
the set {1, 2, ...n}. The point is of course that the j-th hit can be identical to 
the i-th hit, so this is the well-known problem of an ordered set with repetitions. 
The number of different such fc-fold hit events is then given by the formula: 

Z k , n = n k (20) 

where we introduced the symbol Zk, n to stand for the number of ordered sets of 
fc hits out of, each time, n possibilities. Let us remind, for completeness, that, if 
we were not to allow repetitions (that is, we require the fc hits to be different), 
that the number of possibilities is the number of permutations: 

p i - (21) 

We now want to find out, in the case of ordered sets with repetitions, what 
is the number of different fc-fold hit events, in which there are exactly I different 
values occurring in the series, where of course I < fc and / < n. Indeed, because 
of the possibility of repetitions, some of the fc different drawings can be identical, 
and we want to know how many different cells are finally hit by the fc hits. We 
will call this number: Z l k n . We didn't find any explicit reference to this number 
in the literature, but one can easily find a recursive definition for it. First, we 
consider special cases: 

• The case I = fc. We want all fc hit cells to be different. This is just drawing 
without repetition: 

Z k kn = Pk = 7 "L, (22) 
fc ' n fe (n-fc)! y ' 

• The case 1 = 1. All fc hits are identical. Clearly, there are exactly n 
different ways to do so: 

Zln = n (23) 

• The recurrence step. In order to compute Z l k n , we consider that we have 
already (fc — 1) hits, and we add the last hit. There are two possibilities: 
or there are already / different elements hit in the (fc — 1) first hits, in 
which case the fc-th hit must be one of these I elements ; or there were 
only (I — 1) different elements hit in the first (fc — 1) hits, and hence the 
last hit must be a different hit, drawn from the (n — l + 1) remaining un-hit 
cells. As such, we obtain: 



J k,n 



Z{- hn .l + Z\-\ n .(n - I + 1) (24) 
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The two special cases and the recurrence relation define the function Z\ n com- 
pletely: indeed, at each recurrence step, k is diminished by 1, while I remains, 
or is diminished by 1. So sooner or later, k will reach the value of I, in which 
case the first special case is applicable, or I will reach the value of 1, in which 
case the second special case is applicable. 
One should obtain that: 

k 

J2zl n = Z k , n =n k (25) 
i=i 

Indeed, when summing over all possibilities of exactly I different cells, we 
obtain of course all the possibilities of drawing k hits out of n possible cells. We 
didn't succeed in obtaining an entirely closed form for the solution, but we did 
find solutions in closed form for small values of 1 0: 

zln = (-2 + 2^)(-l + B )» (2g) 
(27 - 272 1+fc + 173 fc ) (-2 + n) (-1 + n) n 

^k,n ~ (/'J 

4 _ (-384 + 9 2 8 +' £ - 2176 3 fc + 555 4 fe )(7i-3)(n-2)(?i-l)n 
Z k , n - ^ (28) 



A. 2 Counting indistinguishable hits. 

Consider, for completeness^, the case where the hits are indistinguishable in 
principle. In that case, the order of the hits in a fc-hit event is unimportant. 
We can introduce a similar quantity to Z kn , namely, K k n , with the meaning 
of I different cells being hit in a fc-hit event out of n possible different cells, 
but this time where the order is unimportant (though the multiplicity is). The 
total number of different (unordered) multisets with repetition with cardinality 
k drawn out of n possible values is known to be: 

(;>-( b+ **- 1 )-( b ;!7 1 ) < 29 > 

1 These relations are found iteratively by solving the single- variable recursion relation in k 
for respectively I = 2, I = 3, I = 4 ; however, in order to write out this recursion relation for I, 
we need already the explicit expressions for i— 1. For instance, for 1 = 2, the recursion relation 
is found to be: Z\ n = n(n — 1) + 2Z%_ 1 , which is a pure recursion relation over k. We 
can solve it by standard difference equation techniques. Once we know Z? , we can write a 
recursion relation for Z| n = (n — 2)Z'^_ 1 n +3Z^,_ 1 n in which the underlined part is replaced 
by its (now known) explicit expression, so that only an explicit recursion in k remains. In 
all generality, we have: Z l k n = (n — I + n + IZ\,_ X n , again, with the underlined part 

replaced by its explicit closed-form expression. 

2 In usual wire chambers, the quantum-mechanical indistinguishability of multi-particle 
states doesn't play a role as the phase-space resolution given by the "pixelisation" in space and 
time of the detector is several orders of magnitude more coarse than the quantum-mechanical 
"pixelisation". Nevertheless, in the case where the bosonic character of the particle states 
plays a role, we only have to replace Z by K. 
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Here, ( ^ J = is the binomial coefficient. 

If we consider k unordered hits without repetition then we simply have the 
number of combinations, given by the binomial coefficient. From this, we can 
deduce that: 

Kt,n = C k n (30) 

We can also deduce that: 

i2 K L = ( n k ) (31) 

1=1 

In order to derive a general expression for K k , we consider the following 
ordering of the k hits: the first I cells are the cells which have distinct hits, and 
the last k — I cells are the cells which have already been hit. The first I cells (all 
different) can be chosen in C n different ways out of the n existing cells. The 
last k — I cells have to be all the unordered drawings with repetition out of the 

I cells that have already been chosen in the first part: there are ( , ^ , ) ways 

of doing this. So we conclude that: 

Kk,n = (32) 

and we have an explicit, closed- form expression for the unordered multisets of 
cardinality k, with exactly I distinct elements, drawn out of n possibilities. 



B Dead times. 

The classical models of paralyzable and non-paralysable detectors (see [1] for 
instance) are defined as follows: a non-paralyzable detector is "blind" for a fixed 
time r after an event is registered (no matter what happens during this time), 
and is sensitive again after this time. The relationship between the true rate n 
and the observed rate m is then given by: 

(33) 



1 — rnr 



If the detector is "far from saturation" (relatively few events are lost), then we 
can write: 

m ~ n(l — rn) (34) 

In the case of a paralyzable detector, the detector is "blind" for at least a time 
r after an event is registered, but each time there is a new (uncounted) event 
during the "blind period", the detector remains blind for a time r after this 
event. In other words, the detector becomes ready again only if there is a 
time of "silence" after an event of at least r. A detectable event must hence 
be preceded by a silence of at least r. We summarise the derivation given in 
PP, because it will be instructive for our own developments in this text. The 
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probability, given an event at time t = 0, to have the next event happening 
between time t and time t + dt, equals, for a Poisson stream: 

I\(t)dt — nexp(—nt)dt (35) 

The probability to have an event occurring after a "silence" of at least time r 
is then given by 

/>oo 

(36) 



/■OO 

/ Ii(t)dt = cxp(— nr) 

Jt = T 



This gives then: 

m = nexp(— nr) (37) 
If the detector is "far from saturation", then we can write also: 

m ~ n(l — Tin) (38) 

so, far from saturation, there is no difference in behavior between a paralyzable 
and a non-paralyzable detector. 
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